Search for: All records

Abstract Exotic annual grass invasions in water‐limited systems cause degradation of native plant and animal communities and increased fire risk. The life history of invasive annual grasses allows for high sensitivity to interannual variability in weather. Current distribution and abundance models derived from remote sensing, however, provide only a coarse understanding of how species respond to weather, making it difficult to anticipate how climate change will affect vulnerability to invasion. Here, we derived germination covariates (rate sums) from mechanistic germination and soil microclimate models to quantify the favorability of soil microclimate for cheatgrass (Bromus tectorumL.) establishment and growth across 30 years at 2662 sites across the sagebrush steppe system in the western United States. Our approach, using four bioclimatic covariates alone, predicted cheatgrass distribution with accuracy comparable to previous models fit using many years of remotely‐sensed imagery. Accuracy metrics from our out‐of‐sample testing dataset indicate that our model predicted distribution well (72% overall accuracy) but explained patterns of abundance poorly (R² = 0.22). Climatic suitability for cheatgrass presence depended on both spatial (mean) and temporal (annual anomaly) variation of fall and spring rate sums. Sites that on average have warm and wet fall soils and warm and wet spring soils (high rate sums during these periods) were predicted to have a high abundance of cheatgrass. Interannual variation in fall soil conditions had a greater impact on cheatgrass presence and abundance than spring conditions. Our model predicts that climate change has already affected cheatgrass distribution with suitable microclimatic conditions expanding 10%–17% from 1989 to 2019 across all aspects at low‐ to mid‐elevation sites, while high‐ elevation sites (>2100 m) remain unfavorable for cheatgrass due to cold spring and fall soils. 

Abstract: One strand of modern coexistence theory (MCT) partitions invader growth rates (IGR) to quantify how different mechanisms contribute to species coexistence, highlighting fluctuation‐dependent mechanisms. A general conclusion from the classical analytic MCT theory is that coexistence mechanisms relying on temporal variation (such as the temporal storage effect) are generally less effective at promoting coexistence than mechanisms relying on spatial or spatiotemporal variation (primarily growth‐density covariance). However, the analytic theory assumes continuous population density, and IGRs are calculated for infinitesimally rare invaders that have infinite time to find their preferred habitat and regrow, without ever experiencing intraspecific competition. Here we ask if the disparity between spatial and temporal mechanisms persists when individuals are, instead, discrete and occupy finite amounts of space. We present a simulation‐based approach to quantifying IGRs in this situation, building on our previous approach for spatially non‐varying habitats. As expected, we found that spatial mechanisms are weakened; unexpectedly, the contribution to IGR from growth‐density covariance could even become negative, opposing coexistence. We also found shifts in which demographic parameters had the largest effect on the strength of spatial coexistence mechanisms. Our substantive conclusions are statements about one model, across parameter ranges that we subjectively considered realistic. Using the methods developed here, effects of individual discreteness should be explored theoretically across a broader range of conditions, and in models parameterized from empirical data on real communities. 

Variability of the terrestrial global carbon sink is largely determined by the response of dryland productivity to annual precipitation. Despite extensive disturbance in drylands, how disturbance alters productivity-precipitation relationships remains poorly understood. Using remote-sensing to pair more than 5600 km of natural gas pipeline corridors with neighboring undisturbed areas in North American drylands, we found that disturbance reduced average annual production 6 to 29% and caused up to a fivefold increase in the sensitivity of net primary productivity (NPP) to interannual variation in precipitation. Disturbance impacts were larger and longer-lasting at locations with higher precipitation (>450 mm mean annual precipitation). Disturbance effects on NPP dynamics were mostly explained by shifts from woody to herbaceous vegetation. Severe disturbance will amplify effects of increasing precipitation variability on NPP in drylands. 

Ecosystems are experiencing changing global patterns of mean annual precipitation (MAP) and enrichment with multiple nutrients that potentially colimit plant biomass production. In grasslands, mean aboveground plant biomass is closely related to MAP, but how this relationship changes after enrichment with multiple nutrients remains unclear. We hypothesized the global biomass–MAP relationship becomes steeper with an increasing number of added nutrients, with increases in steepness corresponding to the form of interaction among added nutrients and with increased mediation by changes in plant community diversity. We measured aboveground plant biomass production and species diversity in 71 grasslands on six continents representing the global span of grassland MAP, diversity, management, and soils. We fertilized all sites with nitrogen, phosphorus, and potassium with micronutrients in all combinations to identify which nutrients limited biomass at each site. As hypothesized, fertilizing with one, two, or three nutrients progressively steepened the global biomass–MAP relationship. The magnitude of the increase in steepness corresponded to whether sites were not limited by nitrogen or phosphorus, were limited by either one, or were colimited by both in additive, or synergistic forms. Unexpectedly, we found only weak evidence for mediation of biomass–MAP relationships by plant community diversity because relationships of species richness, evenness, and beta diversity to MAP and to biomass were weak or opposing. Site-level properties including baseline biomass production, soils, and management explained little variation in biomass–MAP relationships. These findings reveal multiple nutrient colimitation as a defining feature of the global grassland biomass–MAP relationship. 

Abstract Matrix population models are frequently built and used by ecologists to analyse demography and elucidate the processes driving population growth or decline. Life Table Response Experiments (LTREs) are comparative analyses that decompose the realized difference or variance in population growth rate () into contributions from the differences or variances in the vital rates (i.e. the matrix elements). Since their introduction, LTREs have been based on approximations and have not included biologically relevant interaction terms.We used the functional analysis of variance framework to derive an exact LTRE method, which calculates the exact response of to the difference or variance in a given vital rate, for all interactions among vital rates—including higher‐order interactions neglected by the classical methods. We used the publicly available COMADRE and COMPADRE databases to perform a meta‐analysis comparing the results of exact and classical LTRE methods. We analysed 186 and 1487 LTREs for animal and plant matrix population models, respectively.We found that the classical methods often had small errors, but that very high errors were possible. Overall error was related to the difference or variance in the matrices being analysed, consistent with the Taylor series basis of the classical method. Neglected interaction terms accounted for most of the errors in fixed design LTRE, highlighting the importance of two‐way interaction terms. For random design LTRE, errors in the contribution terms present in both classical and exact methods were comparable to errors due to neglected interaction terms. In most examples we analysed, evaluating exact contributions up to three‐way interaction terms was sufficient for interpreting 90% or more of the difference or variance in .Relative error, previously used to evaluate the accuracy of classical LTREs, is not a reliable metric of how closely the classical and exact methods agree. Error compensation between estimated contribution terms and neglected contribution terms can lead to low relative error despite faulty biological interpretation. Trade‐offs or negative covariances among matrix elements can lead to high relative error despite accurate biological interpretation. Exact LTRE provides reliable and accurate biological interpretation, and the R packageexactLTREmakes the exact method accessible to ecologists.